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Abstract 

We present the general rules for double- Reggeon production of objects with different 
spins and parities. The existing experimental information on resonance production in 
these central exclusive diffr active processes is discussed. The absorptive corrections are 
calculated and found to depend strongly on the quantum numbers of the produced states. 
The central exclusive diffractive production of 0+ and 0~ Higgs bosons is studied as an 
illustrative topical example of the use of the general rules. The signal for diffractive O'^ 
and 0^ Higgs production at the LHC is evaluated using, as an example, the minimal 
supersymmetric model, with large tan /?. 



1 Introduction 

It is always a challenge to measure the quantum numbers of new states, particularly their 
spin and parity. We may measure the specific characteristics of given decay channels or angular 
correlations of the accompanying particles in the production process, especially the correlations 
between the outgoing protons in the central exclusive production process, pp — > p+h+p, shown 
in Fig. 1(a). The advantage of the latter approach is that it offers the possibility, not only 
to separate different states by accurately measuring the missing mass, but also to distinguish 
between scalar and pseudoscalar new heavy objects, which is difficult from studying their decay 
products. In this paper we begin by studying the general implications of applying Reggeon 
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Figure 1: (a) The central production of a state h by double-Reggeon exchange, (b) The double- 
Pomeron exchange contribution io pp p + h + which dominates at high energies, where 
the + signs are used to indicate the presence of Pomeron-induced rapidity gaps. 

techniques to describe such exclusive processes. At very high energies, and in the central region 
[xp — 0), the double- Pomeron process. Fig 1(b), should give the dominant contribution. 

The theory of double-Reggeon (and multi-Reggeon) exchanges was developed long ago [1]. 
However the revival of interest in these processes is related to the new effects observed in the 
central production of resonances by the WA102 experiment [2] in the reaction pp — >■ pX^p, and 
to the proposal to look for the Higgs boson and other new particles in double- Pomeron-exchange 
processes, see, for example, [3]-[7]. Indeed it will be one of the main challenges of the LHC to 
identify the nature (including the spin-parity) of newly-discovered heavy objects. It appears 
to be very hard to find a spin-parity analyser using conventional approaches. 

Models for double-Pomeron-exchange production of hadrons with different quantum num- 
bers have been developed in recent years [8]-[13]. However in some papers [9, 10] the formulas 
of Reggeon theory were not fully consistent, while some results of the others follow simply from 
general rules of the Reggeon approach. In Section 2 we first consider these rules and compare 
them with experimental observations [2], and with the results of the phenomenological anal- 
ysis performed in Ref. [12]. Also the dynamics for the Pomeron-Pomeron-particle vertices is 
discussed. 

In Section 3 we illustrate how the general behaviour is distorted by the dynamics of the 
process, using /i(O^) exclusive diffractive production as an example. Apart from subsection 3.1, 
where we discuss the uncertainties in the predicted cross sections, this section neglects the 
absorptive or unitary corrections. However at high energies these corrections are important 
(see, for example, Refs. [14]-[18]). They correspond to the diagrams of Fig. 2, and can be 
calculated using the Reggeon diagram technique [19]. It will be shown in Section 4 that the 
inclusion of these diagrams leads not only to a decrease of the cross sections of the double- 
Pomeron processes, but also to significant modifications of the angular correlations between the 
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Figure 2: Unitarity or rescattering corrections to the exclusive diffractive process pp 



p+h+p. 



outgoing (forward) protons. Moreover, the magnitude of these absorptive corrections depends 
on the quantum numbers of the produced state h. 

In Section 4 we illustrate the results using the important topical example of the double- 
Pomeron production of heavy bosons. In particular we compare the production of a Standard 
Model Higgs boson with spin-parity = 0~^, with that for a 0~ Higgs^ which appears in various 
extensions of the Standard Model, in particular in a supersymmetric extension. In Section 5, 
we consider the consequences of this approach to investigations of the Higgs sector at the LHC. 
For illustration we evaluate the exclusive cross sections using the minimal supersymmetric 
model^ (MSSM) with large tan/5; a domain in which, for m(O^) ^ 200 GeV, the conventional 
searches at the LHC will face difficulties to discriminate between the different Higgs states and 
to determine their masses. This is especially true in the so-called "intense coupling limit", 
rrih ~ rriA ~ rriH ~ (9(100 GeV) [21]. As a specific example we calculate the event rates for the 
exclusive central production of mass 115 GeV 0^ bosons at the LHC. 



2 Exclusive diffractive production: general rules 

Here we study the general rules for the amplitudes of the exclusive diffractive process 

1 + 2 3 + /1 + 4, (1) 

shown in Fig. 1, where 1, ... ,4 are hadrons, and where the centrally produced particle h has 
spin and parity J^. We show that the production process has characteristic features, that 
depend on the value of , which follow from general principles. 

^For convenience of presentation we will denote this particle h{0^), rather than the conventional notation 
A(O-). 

^For a recent review see, for example, Ref. [20]. 
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To begin, we assume that all the particles are spinless. Then, at high energies and small 
momentum transfer, 

Sl = {P3 + Phf > So, S2 = (P4 + PhY > So 

^1 - -P3± ~ So, t2 ~ -pI± < So, 



T^2\si,S2,P3±,M = E ^13(^1)^24(^2) - - vMti))ri{akit2))g1k{h,t2,cl>), 

ik ^'^0^ 



the amplitude can be written in the form [1] 

.SqJ \So. 

(3) 

where sq = 1 GeV^, 4> is the angle between the transverse momenta ps± and p4± of the outgoing 
protons and 

is the signature factor for Regge pole i with trajectory ai{t) and signature (jj = ±1. The vertex 
factors ^'13(^1) and 5^24(^2) describe the 13 — ai and 24 — ak couplings respectively, while g^j^ 
describes the transition ajCCfc — > h. Note that g^f. depends, in general, on all the scalars that 
can be formed from the vectors which enter the vertex. Moreover, in the case of Rcggcons, the 
longitudinal and transverse components of their momenta act as two different vectors [19]. In 
our case, where the mass of the boson h is fixed, it is enough to keep the transverse momenta 
'P3± and p4±, and the unit vector hq in the direction of the colliding hadrons. Unlike gis and 
5124, the function g^i^ may be complex. In the Regge domain (2), 

S1S2 = s(m^+p^J, (5) 

where s = {pi +^2)^- 

When spin is included, the process is described by helicity amplitudes, each of which has a 
double- Regge behaviour as in (3) [22]. 

T^S.' {^i,S2,h,t2,<f>) = J29x,Xs{ti)gx,xAt2)i-) - vMh))r)iak{t2))gtHti,h,<f>) 

ik VSo/ VSo/ 

(6) 

The relations between the vertex couplings for different helicities, due to conservation of parity 
and other quantum numbers, are the same as for 2 — > 2 reaction [23]. For example, 

g,,,,{t) = {-l)^^-'^^^g^x,-xs{t), (7) 

with 

^l^VM-lf'-^'Pi'Ji: (8) 

where 5*1 {S3) and r/i (773) are the spin and parity of the particle 1 (3) respectively, and Pj, CTj 
are the parity and signature of the Reggeon i. The vertices behave as [23, 24], 

gxMt)-^{-t)^''-''^^', ast^O. (9) 
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Note that relation (7) depends only on the product PjUj and thus the model-independent spin 
structure of the vertices 5'AiA3(fl'A2A4) is the same for all Reggeons with the same product PiCTj. 

Below we will be interested mainly in the spin structure of the central vertex gl^{t\, ^2,0)- 
It can be written in the form [22] 

oo 

9ik'= E e^'^^S^U' withmi+m2 = A,, (10) 

m2=—oo 

where m2 has the meaning of the projection of the angular momentum jk of Reggeon k (ana- 
lytically continued from all angular momenta jk in the , t2-channels) . Now invariance under 
parity leads to the relation [22] 

^(_i)\ht^-\h (II) 

im\m2 \ J /— mi— m2' v / 

where 

= rih{-lf'^Pi<yiPkak. (12) 

Thus the spin structure of the central vertex depends only on the product of the naturalities 
(that is the parities and signatures) of particle h and the exchanged Reggeons^. The behaviour 
of '^^im2 small ^1,^2 is given by the formula [22] 

ltm2 ~ With mi + m2 = A,. (13) 

Note that all values of m2 (mi) enter (10), but, due to (13), for small ti,t2 it is enough to 
consider the lowest values of m2 (mi) consistent with (11). 

It is often convenient to write the spin structure of the amplitudes in terms of the char- 
acteristic 3-vectors of the problem. Such a representation is closely related to the hehcity 
amplitudes discussed above [24], but the formulas become more transparent. In this case the 
central vertex g^f. is written as a scalar (or pseudoscalar) function (depending on the product 
r]h{—lY''(TiPi(TkPk) of the vectors P3±,P4± and the spin vectors (tensors) of particle h. 

Let us consider particular examples for the spin-parity of /i, in each case taking 
(TiPiUkPk = +1 as for double- Pomeron exchange. 

(a) jP{h) = 0+ 

For a scalar particle /i, the vertex coupling is simply 

9ik = fo+ ipl± > pI± > P3± ) , (14) 

where /o+ is a function of the scalar variables which can be formed from the transverse 
momenta p3± and p4± of the outgoing protons. When or p4± 0, this function in 
general tends to some constant /(0,0,0). In order to obtain further information on the 
structure of this function, extra dynamical input is needed (sec Section 3). 

^It can be shown that similar formulae are also valid for photon-photon fusion processes. 
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(b) J^(/i) = 0- 

For the central production of a pseudoscalar particle, the vertex factor takes the form 

fife = /o-(P3±>P4±>P3±-P4±) £ikl{P3±)i{P4±)k{no)l, (15) 

where no is the unit vector in the direction of the colliding hadrons (in the c.m.s.) . In this 
case, according to (11), all amplitudes with mi,m2 = are zero, and so |mi| = |m2| = 1 
are the leading terms. According to (13) this corresponds to helicity amplitudes which 
are proportional to (— ti)2(— ^2)2 for small t. Thus the cross section behaves as ~ |ti||t2|- 
Also the angular distribution contains a factor sin^ 0, which is evident from either (15) 
or (10) and (11). Again the function /q- is not predicted from general principles. 

Note that the characteristic sin^ (f) dependence of the angular distribution, and the t- 
behaviour at small t, which are observed by the WA102 Collaboration for 77, and rj' 
production [2], are direct consequences of the general properties of the double-Rcgge- 
exchange amplitudes. This behaviour is valid not only for doublc-Pomcron exchange, 
but also for Pf,ff,pp,A2A2,uu, . . . exchanges. Under the interchange of the Reggeons, 
i k, 

{no)i -{no)i (16) 

and hence if the Reggeons are the same, i = k, the function /q- should be symmetric 
under the interchange 3 <-> 4. 

(c) r{h) = 1+ 

For the production of an axial vector state, both — and = ±1 components are 
present, and the vertex factor can be written as 

9ik = fi+ £iki{P3±)i{P4±)kei + (/iV {P3±)i + 11+ (P4±)i) eiki{no)kei, (17) 

where the /* are functions of the scalar variables pl±,p'i± and p3± ■ p4±, and where e 
is the polarization vector of the 1"*" meson. If both of the exchanged Reggeons are the 
same {i = k), then the function g^^ is symmetric under the interchange 3 4. As a 
consequence, for small Pi±, 

/^-(pL-pL), = (18) 

The form of the vertex factors g^^. for the central production of states h of higher spin can 
readily be constructed using equations (10)-(12). 

It is interesting to note that the structure of the vertex factors given in (14)-(18) coincides, 
for small Pi±, with that found using a non-conserved vector current model [12], which gives a 
good description of the experimental data of the WA102 Collaboration [2]. However, from the 
discussion above, it is clear that these results simply follow from the general rules of Reggeon 
theory. They are not connected with a particular vector current model of the Pomeron, but 
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rather follow from the fact that the product of the parity and signature of the Pomeron is 
+1. Moreover the Pomeron has positive signature and corresponds to the analytic continuation 
from angular momenta = 2+, 4+, . . . in the i-channel. The same results would be obtained 
if a tensor current model were used. 

On the other hand, the detailed structure of the amplitudes fm{P3±: P4±j P3±'P4^±) (with m = 
0^, 0^, 1^, 2+, . . .) depends on dynamics and cannot be predicted from the general principles of 
Regge theory. For example, if the heavy state h is strongly coupled to gluons and is produced 
perturbatively via the diagram shown in Fig. 3(a), then in the forward direction {p3±,p4± -C Q±) 
the vertex factor /o+ does not depend on ps or p4. Moreover, as was shown in [25, 26, 7], 
there exists a — 0, parity-even, selection rule, for production by gluon-gluon fusion where 
each of these active gluons comes from colour-singlet digluon ^-channel exchange, see Fig. 3. 
As a consequence the production of the negative-parity h{0~) state is strongly suppressed 
in comparison with the production of the h{(}~^) state. Similarly it follows that the central 
diffractive exclusive production of 2++ states is also suppressed in some topical cases; for 
example, 2^^ states formed from heavy quark pairs (in the non-relativistic model) [26, 27] or 
2++ 'gravitons' in models with extra dimensions in which their coupling to gluons has a point- 
like nature (with no derivatives) so they arc not produced via the = two-gluon state [7]. 
Also these processes can provide a unique opportunity to determine the quantum numbers of 
pair-produced new strongly-interacting objects [7]. For example, comparatively light gluinos 
and squarks can be distinguished by the respective and /3 threshold behaviour, where (3 is 
the particle velocity. 



3 Example: dynamics of /i(O^) Higgs production 

So far we have discussed the structure of the production amplitudes for 

pp^p + h + p, (19) 

where h has a given J^, which follow from general principles. To go further we need to 
consider the dynamics of the process. We study h{0'^) Higgs production as a topical example. 
The general rules imply that the central vertices behave as 

9ik°'^^ ~ constant , , 

9ik ~ (P3± X P4_l) • no ~ |ti|2|t2|2 sm0 

at small t. 



3.1 Amplitudes for h{0^) production 

To see how the dynamics modify this behaviour we have first to describe how the cross sections 
for the exclusive production of /i(O^) Higgs bosons are calculated. We use the formalism of 
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Figure 3: (a) The QCD diagram for double-diffractive exclusive production of a Higgs boson 
h, pp —>■ p + h + p, where the gluons of the hard subprocess gg —>■ h are colour screened by 
the second t-channel gluon. (b) The rescattering or absorptive corrections to pp p + h + p, 
where the shaded region represents the multi-Pomeron exchanges of Fig. 2. 



Ref. [6, 26, 7]. The amplitudes are described by the diagram shown in Fig. 3(a), where the hard 
subprocesses gg — > h{0^) are initiated by gluon-gluon fusion and where the second t-channel 
gluon is needed to screen the colour flow across the rapidity gap intervals. The Born amplitudes 
are of the form 

Tk = N f ^'^^J^ — -/,(xi,x;,Q^,/.^;ti)/,(x2,x'2,Q^,/i^;t2) (21) 

The overall normalization constant can be written in terms of the h ^ gg decay width [7], 
and the gg ^ h vertex factors are 

^"^(0+) = {Q± - P3±) ■ {Q± + P4±) . s 

Vh{o-) = ({Q±-P3±) ^ {Q±+P4±)) -no. 



The /g's are the skewed unintegrated gluon densities of the proton at the hard scale fi, taken 
to be m/j/2, with 

Qs = minlQ^, KQ^ -P3_l) 
Q4 = minjQ^, KQ^ +P4_l) 

Below, we assume factorization of the unintegrated distributions. 



(23) 



fg{x,x',Q\fi';t) = fg{x,x',Q\fi')FN{t), (24) 

where we parameterize the form factor of the proton vertex by the form -F/v(t) = exp(6t) with 
6 = 2 GeV~^. To single log accuracy, we have [28] 

d 



f,ix, x\ Ql /x^) = R^Q^ ( Vr(g., /x) xg{x, Ql) ) , (25) 



8 



where T is the usual Sudakov form factor which ensures that the gluon remains untouched in 
the evolution up to the hard scale so that the rapidity gaps survive. The square root arises 
because the bremsstrahlung survival probability T is only relevant to the hard gluon. Rg is the 
ratio of the skewed x' <^ x integrated distribution to the conventional diagonal density g{x^ Q"^). 
For a; ^ 1 it is completely determined [29]. The apparent infrared divergence of (21) is nullified^ 
for h{0'^) production by the Sudakov factors embodied in the gluon densities fg. However the 
amplitude for /i(0 ) production is much more sensitive to the infrared contribution. Indeed 
let us consider the case of small pi± of the outgoing protons. Then, from (22), we see that 
V/i(o+) ~ Q\, whereas T4(o-) ~ P3±P4± (since the linear contribution in Q± vanishes after the 
angular integration). Thus the (PQ±/Q\^ integration for /i(0+) is replaced by P3±P4:±(PQ±/Q± 
for h{0~), and now the Sudakov suppression is not enough to prevent a significant contribution 
from the Qi ^ 1 GeV^ domain. 

3.2 Uncertainties 

To estimate the uncertainty in the predictions for the h'^iO) exclusive diffractive cross sections 
we first quantify the above uncertainty arising from the infrared region, where the gluon dis- 
tribution is not well known. Fig. 4 shows the dependence of h{0^) and /i(0+) production at 
the LHC, for ruh — 120 GeV and /i = mh/2, using different treatments of the infrared region. 
The continuous and dashed curves are calculated using MRST99 [30] and CTEQ6M [31] par- 
tons respectively with the very low Q gluon frozen at its value at (53,4 = 1.3 GeV. Then we 
integrate down in Q± until (^3^4 are close to Aqcd, where the contribution vanishes due to the 
presence of the T-factor. This will slightly overestimate the cross sections as the gluon density 
decreases with decreasing for a; ~ 0.01. A lower extreme is to remove the contribution below 
Q3,4 = 1.3 GeV entirely. The result for MRST99 partons is shown by the dotted curve in Fig. 4. 
Even with this extreme choice, the 0+ cross section is not changed greatly; it is depleted by 
about 20%. On the other hand, as anticipated, we see for 0~ production, the infrared region is 
much more important and the cut reduces the cross section by a factor of 5. 

Another uncertainty is the choice of factorization scale fi. Note that in comparison with 
previous calculations [7], which were done in the limit of proton transverse momenta, P3,4_L -C 
Q±, now we include the explicit ]9_L-dependence in the Q^-loop integral of (21). Moreover, we 
resum the 'soft' gluon logarithms. In 1/(1 — 2;), in the T-factor.^ So now the T-factor includes 
both the single soft logarithms and the single coUinear logarithms. The only uncertainty is the 
noii-logaritriniic NLO eontiil)ution. This may be modelled by changing the factorization scale, 

**In addition, at LHC energies, the effective anomalous dimension of the gluon gives an extra suppression of 
the contribution from the low Q± domain [4]. 

^To account for the interference and the precise form of the amplitude for soft gluon {q± <C irifi) large-angle 
emission, we explicitly calculate the one-loop virtual correction to the gg ^ h vertex, integrating over the whole 
angular range for qj_ <C mh- We adjust the upper limit of the z-integral so that z < 0.62m;i/(0.62TO;i + q±) in 
the expression for the T-factor (see eq. (10) of [7] with = q±), in order to reproduce the complete one-loop 
result. 
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Figure 4: The 4> dependences of diffractive exclusive h{0'^) and ^(0^) = h{0~) production, 
pp — > p + /i + p, at the LHC, with rrih — 120 GeV. is the angle in the transverse plane 
between the outgoing protons. The curves are for central rapidity and do not include absorptive 
corrections. They correspond to the effective luminosities for the gg — > h{0'^) subprocesses 
integrated over the outgoing proton momenta Pi±, assuming an exp(— 6pfj^) behaviour of the 
unintegrated gluon densities fg, with slope b — 2 GeV~^. 



/X, which fixes the maximal q± of the gluon in the NLO loop correction. As the default wc have 
used II = mh/2; that is the largest g_|_ allowed in the process with total energy rrih. Choosing a 
lower scale /x = would enlarge the cross sections by about 30%. 

Next there is some uncertainty in the gluon distribution itself. To illustrate this, we compare 
predictions obtained using CTEQ6M [31], MRST99 [30] and MRST02 [32] partons. For 0+ 
production at the LHC, with nih — 120 GeV and /x = mh/2, we find that the effective gluon- 
gluon luminosity, before screening, is 



dC 



dydlnM^ 



= (2.2, 1.7, 1.45) X 10-2 (26) 
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respectively. This spread of values arises because the CTEQ gluon is 7% higher, and the 
MRST02 gluon 4% lower, than the default MRST99 gluon, in the relevant kinematic region. 
The sensitivity to the gluon arises because the central exclusive diffractive cross section is 
proportional to the 4th power of the gluon. For 0~ production, the corresponding numbers 
are (4.2, 2.7, 1.7) x 10^^. Up to now, we have discussed the effective gluon-gluon luminosity. 
However, NNLO corrections may occur in the gg ^ h fusion vertex. These give an extra 
uncertainty of ±20%. Note that we have already accounted for the NLO corrections for this 
vertex [7]. 

Finally, we need to consider the uncertainty in the evaluation of the soft rescattering cor- 
rection factor S^, which is the probability that the rapidity gaps survive the soft pp interaction. 
The computation of S'^ is discussed in some detail in Section 4. Here it suffices to say, from the 
analysis [16] of all soft pp data, that a conservative error on the values of is ±50%. 

Combining together all these sources of error we find that the prediction for the 0^ cross 
section is uncertain to a factor of almost 2.5, that is up to almost 2.5, and down almost to 
1/2.5, times the default value.® On the other hand, 0^ production is uncertain by this factor 
just from the first (infrared) source of error, with the remaining errors contributing almost 
another factor of 2.5. Although the rate of h{Q~) production is very sensitive to the infrared 
contribution, and could indicate the presence of a significant non-perturbative contribution, we 
find that the form of the dependence is not. We discuss this point in Section 5. 

Note that the non-local structure of the amplitude leads to an extra angular dependence 
coming from the correlations between Q± and the Pi± in the integral in (21). In fact, expanding 
the gluon propagators gives corrections of the type 

Q±-P3±/Ql, Q±-P4±/Ql, (27) 

which lead to an additional contribution of the form —p3±-p4±/Q\- This reflects the dependence 
of the vertex factors /o± on p3± ■ p4±, see (14) and (15). However, it is evident from Fig. 4 
that this contribTition does not give a large effect. What is more important is the suppression 
of h{0~) production in comparison to that for h{0'^). The h{0~) amplitude is proportional 
to p3± X p4_L, where the dimensions must be compensated by some scale. In perturbative 
QCD this is the scale Q\ arising from the gluon loop in Fig. 3(a). Therefore the h{0^) cross 
section is reduced by a factor {p3±pI±/'2Q'\_ ), that is by a factor of the order of 500 for typical 
- 1/26 ~ 0.25 GeV^ ff Qi ~ 4 GeV^. 

4 Absorptive corrections 

In this section we consider how exclusive double-diffractive production is influenced by the 
absorptive (shadowing) effects, which arise from the multi-Pomeron diagrams of Fig. 2. To 
determine the suppression due to these absorptive corrections, it is convenient to work in 
impact parameter, b, space. 

^For example, we predict the cross section for the exchisive diffractive production of a Standard Model Higgs 
at the LHC, with ruh = 120 GeV, to be in the range 0.9-5.5 fb. 
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4.1 Absorptive effects for h{0'^) production 



The amplitude for the central production of an h{0~^) state, via the double-Pomeron-exchange 
process pp ^ p + h + p, has the form 

T\s, h, h, b) = exp (^~np{s, b^)^ T!^p{s,, s^, h, 62), (28) 

where fip is the contribution of Pomeron exchange to elastic pp scattering in impact parameter 
space 

np(s,b') = ^eM-b'm, (29) 

where cr^ is the Pomeron contribution to the total cross section of the pp interaction, and 

B^lBo + a'pln{s/so) (30) 

is the slope of the elastic pp amplitude. The amplitude Tpp is the Fourier transform, to impact- 
parameter space, of the amplitude T{2^^(si, S2,P3±,P4±) of (1) with i — k — P. That is 

Tj^pisi, S2, K h) = (^)' / d^Ps±dW e'P'^-^'^e-'^^^-^'^ T^p{s,, S2,m,P4±), (31) 

where 6 = 61 + 62 is the Fourier conjugate to g = ps± —p4±. For simplicity, we give the formula 
for a single-channel cikonal, where only intermediate proton states are considered, between the 
Pomeron exchanges in Fig. 2. The extension to the multichannel case is straightforward. In 
the calculations presented here we used the two-channel eikonal model of Ref. [16]. 

Note that if ap{0) — 1 = A > 0, then flp{s,b'^) increases with energy and leads to a 
substantial suppression of cross section at very high energies. Calculations, using the model of 
Ref. [16], show that at Tevatron energies the Born cross section, corresponding to the diagram 
of Fig. 1(b), is suppressed by the multi- Pomeron exchanges of Fig. 2 by a factor of roughly 
0.05. At the LHC the suppression factor' {S^) is 0.026. 

Since the amount of suppression depends on the impact parameter b, it leads to a charac- 
teristic dependence of the factor on the angle between the outgoing protons [18] . This is 

— * 

related to the fact that b is the Fourier conjugate to the vector q = p3± — p^^. If the outgoing 
protons are tagged, then the characteristic peripheral form of the amplitude in 6-space can 
be studied experimentally in double-Pomeron-exchange processes by measuring the dependence 
of the cross section on q. 

We emphasize that the suppression S^, due to absorptive or rescattering corrections, depends 
not only on the particular process, but also on the kinematical cuts which select events in a 
given pi±, (f) domain. Therefore the suppression has to be calculated for each particular 
kinematical configuration. 

""It is interesting to note that the introduction of the Pi± and angular correlations in (21), (22) raise (5^) 
from 0.023 to 0.026. 
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4.2 Comparison of exclusive diffract ive h{0^) Higgs production 

So far wc have discussed absorptive corrections for h{0^) production. Here we compare these 
corrections with those for h{0~) production. For h{0~) production we predict a different q 
behaviour. This originates from (15); the Born double- Pomeron-exchange amplitude for pro- 
cess (1) now contains the kinematical factor (p3j_ x p^jj • no and this, in turn, implies that 

— * — * 

the Fourier transform contains the factor {bi x 62) • ^o- Thus the corresponding amplitude 
,■§2,61,62) tends to zero as bi or 62 ^0. As a result, the suppression arising from 
rescattering is less effective, and the factor S'^ is larger than for h{0'^) production. Also the 
distribution is distorted due to absorption. 

The effect of the absorptive corrections on the angular correlations between the outgoing 
protons was discussed in detail in Ref. [18] for /i(0+) production^. There it was shown that 
the absorptive corrections are largest in the back-to-back configuration where ps± is directed 
against p4±, since in this case both ti ~ —{k± +P3±)^ and t2 — —{k± — P4±)^ can be minimized 

— * 

simultaneously by the same momentum k± transferred in the elastic rescattering, see Fig. 3(b). 
Thus for (f) — 180° the momentum is transferred mainly through the rescattering amplitude. 
The suppression factor S"^ was plotted versus for different choices of ps± and p4± in Ref. [18]. 
It was shown that the diffractive dip (which arises from the maximum cancellation between 
the bare amplitude and rescattering contribution) moves to smaller as the values of Pi± are 
increased. 

Here we calculate S"^ as a function of for h{0~), as well as h{0~^), exclusive diffractive 
production. We integrate over the pi±^ of the outgoing protons assuming an exp{—b{pf^^ +P4±)) 
behaviour of the Pomeron-proton vertices gi^ and 5^24, with 6 = 2 GeV^^. We use the two- 
channel eikonal model of Ref. [16]. For the central vertex we take (^3^ x p4±) ■ hq for h{0~) 
production and a constant for /i(0+) production. The results for the suppression factor S"^ are 
shown in Fig. 5 for h{0'^) production of mass 120 GeV at the LHC energy, ^/s — 14 TeV. As the 
azimuthal angle between p2± and p4± increases, the first diffractive dip, followed by the second 
maximum, are evident in the S'^ curves obtained by integrating over all pi±. The dotted curves 
show the effect of restricting the outgoing protons to the domain pi± < 0.35 GeV. As expected, 
the diffractive dip is pushed to larger angles and is barely reached even for the back-to-back 
configuration, (p = 180°. As we see from Fig. 5 that the survival factor S*^ is about 3-4 times 
larger for h{0^) as compared to h{0^) production. This is a refiection of the more peripheral 
nature of h{0^) production. For the same reason the suppression obtained when integrating 
over the small pi± domain, pi± < 0.35 GeV, is less than when integrating over all pi±, since it 
is more weighted to the larger values of the impact parameter 6. 

Finally in Fig. 6 we show the predictions for the effective luminosity with the absorptive 
effects included. The original sin^ ^ and constant behaviours of h{0~) and /i(0+), respectively, 
are distorted first by the ps± •P4±/Q^ type corrections from the integration over the gluon loop 
in Fig. 4, and then by the absorptive effects given by the suppression factors S"^ shown in Fig. 5. 

^Notc that there is a typographical error in eq. (25) of [18], where the last factor should be simply S"^ rather 
than its second derivative. However the results presented in [18] correspond to the correct definition of F. 
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Figure 5: The suppression factors 5*^ of h{0^) and ^(0^) = h{0~) Higgs production via the 
process pp ^ p + h + p ai the LHC, arising from rescattering effects. The outgoing protons 
are integrated over (i) aU pi± and (ii) pi± < 0.35 GeV (dotted curves). For illustration, the 
continuous curve for /i(0"") = ^(0~) production includes the general sin0 behaviour of the bare 
amplitude. 



5 Consequences for signals in the Higgs sector 

We have studied the central exclusive diffractive production of bosons via the process pp —>■ 
p + h + p, and emphasized that correlations between the outgoing proton momenta reflect the 
spin-parity of h. As a topical example to illustrate these properties we compared h{Q'^) and 
A{0~) = h{0~) Higgs production. In particular, Fig. 6 shows that the dependence on the angle 
between the outgoing proton transverse momenta, ps± and p4±, is different for the natural 
(O"*") and unnatural (0~) parity states of h. The comparison with Fig. 4 shows that absorptive 
effects have significantly distorted the (p distributions and, in fact, have increased the difference 
between the 0+ and 0~ distributions. Thus this distribution provides a unique possibility to 
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Figure 6: The effective luminosities of tlie gg — > A{0 ) = h{0 ) and h{0'^) subprocesses, inte- 
grated over tfie outgoing proton transverse momenta Pi±, in tfie exclusive diffractive processes 
pp ^ p + h + p at the LHC, with ^/s ^ U TeV, using the MRST99 [30] gluon. Absorptive 
(or rescattering) effects are included. The curves are for y = and Higgs masses of 120 GeV. 
Comparison with the continuous curves in Fig. 4 shows the suppression of the event rate, and 
the distortion of the (f) behaviour due to absorptive effects. 



distinguish between 0+ and bosons, which, in the case of inclusive production^, is extremely 
difficult. 

We have seen that the amplitude for the production of unnatural parity (P = (— 1)"^+^) 

®A proposal, similar in spirit to our approach, can be found in Ref. [33]. The idea is to determine the CP- 
parity of a Higgs boson by measuring the azimuthal correlations of the tagging (quark) jets which accompany 
Higgs production via the vector-boson-fusion mechanism. Even if we disregard the possible degradation of the 
characteristic features of the distribution caused by parton showers and the inclusive environment of the jets, 
we note that the method is not applicable in some important regions where the couplings of the Higgs to vector 
bosons are strongly suppressed. Another method to determine the spin-parity of the Higgs, which similarly 
relies on the vector-boson coupling, was discussed in Ref. [34] . 
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states contains a factor (^3^ x p4±) • no. Thus the cross section vanishes as p3± or p4± — > and 
vanishes as sin^ </> as ^ or vr. These properties may be used to suppress the cross section 
for natural parity (P = (—1)'^) production in comparison to that of unnatural parity states. 
In particular, selecting events with P3±,P4± > 0.4 GeV and 20° < < 120° suppresses the 
O"*" yield by about a factor of 10, while only decreasing the 0~ cross section by a factor of 2.3. 
The relative 0~ enhancement may be important as the cross section for the central exclusive 
production of an ^(0~) boson is quite small, and moreover, in many super symmetric scenarios, 
the ^(O"*") (and/or H{0~^)) and A{0~) bosons are close in mass. 

As a specific example we consider Higgs production in the minimal supersymmetric model 
(MSSM) with large tan/3 and ~ 110-130 GeV. In this domain the branching ratios of 
Higgs-like bosons to vector bosons and photon pairs^'^, and the couplings to top quarks, are 
much suppressed [21], and it becomes problematic to perform a complete coverage of the Higgs 
boson sector using the conventional inclusive processes. In particular the problem of resolving 
the signals for different states becomes quite challenging^^. On the contrary the cross section, 
(TcEP, for central exclusive diffractive production in the MSSM is enhanced in comparison with 
that of the SM. The MSSM (for tan (3 = 30) and SM cross sections, o"cep; at the LHC energy, 
are shown in Fig. 7. They have been evaluated using the effective gg luminosities obtained in 
Section 3 with the absorptive corrections calculated in Section 4.2, see Fig. 6. The normalization 
factor in (21) is 



^2 27rr(/i ^ gg) 



where the NLO K factor ~ 1.5 [7] and the number of colours Nc — 3. The widths and properties 
of the Higgs scalar {h, H) and pseudoscalar (A) bosons are calculated using the HDECAY code, 
version 3.0 [35], with all other parameters taken from Table 2 of [35]; also we take IMODEL — 4, 
which means the radiative corrections are included according to Ref. [36]. 

From Fig. 7 we see that the 0+ Higgs bosons should be accessible at the LHC, via the central 
exclusive signal, over a wide mass range up to about 250 GeV in this scenario. The enhancement 
of the MSSM signals is clearly apparent, except near rrih — 127 GeV. For example, for the 
production of a Higgs boson of mass 115 GeV for tan/5 = 30 (or 50) we have 

Br(/i ^ bb) crgp^ ~ 20 (70) fb, (33) 

about 10 (40) times larger than (Jcep in the SM. 

For the same parameters, for A{0~) production we obtain 

Bt{A bb) ac^P ~ 0.2 (0.5) fb. (34) 



^°The branching ratios of h,H,A-^ 77 are less than, or of the order of, 10~^-10~^, which is much smaller 
than in the SM. 

"'^'^The separation of h and H bosons may be especially difficult for inclusive signals, where the mass resolution 
is usually Am ^ 10 GeV, except in the 77 and probably /x/k modes. However, with forward proton taggers, the 
exclusive signal has the added bonus that a mass resolution of Am 1 GeV may be obtained [38] . 
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Figure 7: The cross sections times the bb branching ratio predicted, for the central exclusive 
production of h{0+), H{0+) and ^(0") MSSM Higgs bosons (for tan/? = 30), at the LHC, 
compared with the SM result. We use the MRST99 [30] gluon distribution. The vertical hne 
separates the mass regime of light /i(0+) and heavy -f^(0+) bosons. 



However we have emphasized the infrared sensitivity of the rate of A{0^) production. It is 
possible that a non-perturbativc contribution, coming from low values of in (21), may 
enhance the cross section by a factor of 3 or more. Nevertheless it would be extremely hard to 
observe the ^(0^) boson under the h{0^) signal, when the masses are close. A typical mass 
difference is uia — ruh ^ 2.9 (1.4) GeV for tan (3 — 30 (50), if ttia — 115 GeV. The situation 
for the observation of the ^(0~) boson is even worse due to the comparatively large expected 
widths of the Higgs bosons. For instance, if uia ~ 90-130 GeV and tan (3 — 30, then the widths 
of the Higgs bosons become of order 2 GcV [20] . On the other hand, proton taggers, with a 
very accurate missing mass resolution of AM ~ 1 GcV, offer the attractive possibility, not only 
to separate the h and H bosons, but also to provide a direct measurement of the widths of the 
h (for m/j < 120 GcV) and the H bosons (if > 130 GeV). Also we note that by comparing 
the cross sections of (33) and (34), we see that if a new heavy object were observed in inelastic 
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events, but not in exclusive central production, it would indicate that it had unnatural parity. 

Although the rate of ^(0~) Higgs production is sensitive to contributions from the infrared 
region, we do not expect a significant change in the qualitative behaviour of the 0~ production 
amphtude. The reasons are as follows. First, the vanishing of the amplitude as pi± — > and/or 

— > 0, TT follows from the general form (15) of the vertex in Regge theory. Second, as a rule the 
extra (f) dependence caused by the ■ p4± term is weak, see Fig. 4 for example. Third, in the 
very extreme case where we use GRV94 partons [37] (which enable us to take a very low infrared 
cut-off Qq = 0.4 GeV^, but which are known to overestimate significantly the low x gluon), 
the 0^ Higgs cross section is enhanced, relative to that obtained using MRST99 partons with 
Qo = 1.3 GeV, by about a factor of 4, but the (j) and Pi± dependences are essentially unaltered. 

Returning to /i(0'^) production, we see, for the example of (33), that already for an LHC 
luminosity C = 30 fb^^, about 600 (2000) bosons are produced. If the experimental cuts 
and efficiencies quoted in Ref. [38] are imposed, then the signal is depleted by about a factor 
of 6. This leaves about 100 (400) observable events, with an unaltered background of about 
3 events [38] in a AM = 1 GeV missing mass bin; which gives an incredible significance for a 
Higgs signal! 
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